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TITLE OF THE INVENTION 

RADIUS-IN-IMAGE DEPENDENT DETECTOR ROW FILTERING FOR WINDMILL 

ARTIFACT REDUCTION 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 

[0001] Multi-slice x-ray computer tomography (CT) systems were developed and introduced 
into the medical market around 1998. The number of slices generally ranges from 2 to 40 
currently, and is expected to increase to 64 or even 256. (See Y. Saito, H. Aradate, H. 
Miyazaki, K. Igarashi, and H. Ide, "Development of a Large Area 2-dimensional Detector for 
Real-Time 3-dimensional CT (4D-CT)," Radiology vol. 217(P), 405 (2000); Y. Saito, H. 
Aradate, H. Miyazaki, K. Igarashi, and H. Ide, "Large Area Two-Dimensional Detector 
System for Real-Time Three-Dimensional CT (4D-CT)," Proc. ofSPIEMed. Imag. Conf., 
4320, 775-782 (2001); and http://www3.toshiba.co.jp/medical/4d-ct/, the contents of each of 
which are herein incorporated by reference). One cone-beam image reconstruction algorithm 
with helical scanning for CT systems uses a generalized weighted version of Feldkamp 
reconstruction algorithm, which introduces a flexible focus orbit and applies a weighting 
function to the Feldkamp algorithm. (See L. A. Feldkamp, L. C. Davis, and J. W. Kress, 
"Practical Cone-Beam Algorithm," J. Opt. Soc. Am. A, 6, 612-19 (1984); H. Aradate and K. 
Nambu, "Computed Tomography Apparatus," Japanese Patent No. 2,825,352; L. G. Zeng 
and G. T. Gullberg, "Short-scan Cone Beam Algorithm for Circular and Noncircular Detector 
Orbit," Proc. ofSPIEMed. Imag. Conf, 1233, 453-463 (1990); H. Kudo and T. Saito, "Three- 
Dimensional Helical-Scan Computed Tomography Using Cone-Beam Projections," J. 
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Electron. Information Commun. Soc. Japan, J74-D-II, 1108-1114 (1991); G. Wang, T. H. 
Lin, P. C. Cheng, D. M. Shinozaki, "A General Cone-Beam Reconstruction Algorithm," 
IEEE Trans. Med. Imaging, 12, 486-496 (1993); K. Taguchi, "X-ray Computed Tomography 
Apparatus," U.S. Patent No. 5,825,842 (Filed in 1995); K. Wiesent, K. Barth, N. Novab, et 
al., "Enhanced 3-D -Reconstruction Algorithm for C-arm Systems Suitable for Interventional 
Procedures," IEEE Trans. Med. Imaging, 19, 391-403 (2000); M. D. Silver, K. Taguchi, and 
K. S. Han, "Field-of-View Dependent Helical Pitch in Multi-Slice CT," Proc. ofSPIE Med. 
Imag. Conf., 4320, 839-850 (2001); M. D. Silver, K. Taguchi, and I. A. Hein, "A Simple 
Algorithm for Increased Helical Pitch in Cone-Beam CT," The Sixth International Meeting 
on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, 70- 
73 (2001), the contents of each are herein incorporated by reference). 
[0002] Other reconstruction algorithms include quasi-cone-beam algorithms such as 
advanced single-slice reb inning (ASSR) and adaptive multiple plane reconstruction (AMPR), 
which reconstruct multiple slices that are not aligned to a single axis. The axes of the slices 
fit a helical orbit and may be z-filtered to obtain perpendicular slices. (See Y. Machida, 
"Computed Tomography System," Japanese Patent Disclosure (Kokai) 08-187240; M. 
Kachelriess, S. Schaller, W. A. Kalender, "Advanced Single-slice Rebinning in Cone-Beam 
Spiral CT," Medical Physics vol. 27, pp. 754-772 (2000); S. Schaller, K. Stierstorfer, H. 
Bruder, M. Kachelriess, and T. Flohr, "Novel Approximate Approach for High-Quality 
Image Reconstruction in Helical Cone Beam CT at Arbitrary Pitch," Proc ofSPIE Vol. 4322, 
pp. 113-127 (2001); T. Flohr, K. Stierstorfer, H. Bruder, J. Simon, A. Polacin, and S. 
Schaller, "Image Reconstruction and Image Quality Evaluation for a 16-slice CT Scanner," 
Medical Physics vol. 30, pp. 832-845 (2003), the contents of all of which are herein 
incorporated by reference). 
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[0003] There are some other algorithms which have the same aliasing problem by using 
cone-to-parallel fan-beam rebinning (See, e.g., H. Tuy, "3D Image Reconstruction for Helical 
Partial Cone-beam Scanners," proc of Fully 3D 1999, pp. 7-10; H. Turbell, et al., "Three- 
dimensional Image Reconstruction in Circular and Helical Computed Tomography," 
Licentiate thesis No. 760, Linkoping Univ, ISBN 91-7219-463-4, 1999; H. Turbell, et al., 
"An Improved Pl-method for Reconstruction from Helical Cone-beam Projections," Conf 
record of IEEE MIC 1999; R. Manzke et al, "Extended Cardiac Reconstruction (ECR): A 
Helical Cardiac Cone-beam Reconstruction Method," proc of Fully 3D 2003, Mo-PM2-4., the 
contents of all of which are herein incorporated by reference), and fan-beam algorithms (K. 
Taguchi and H. Anno, "High Temporal Resolution for Multi-slice Helical CT," Medical 
Physics vol. 27, May 2000, the contents of which are herein incorporated by reference). The 
present invention is not limited by the choice of reconstruction scheme. 

DESCRIPTION OF THE RELATED ART 

[0004] One common problem with cone beam or quasi-cone-beam algorithms is insufficient 
sampling intervals in the z-axis (detector row) direction. These algorithms violate the 
Nyquist theorem (which requires two samples within one detector cell aperture) and causes 
aliasing artifacts in reconstructed images due to high frequency components. (See M. Silver, 
K. Taguchi, I. Hein, B. Chiang, M. Kazama, I. Mori, "Windmill Artifact in Multi-Slice 
Helical CT," Proc of SPIE Vol. 5032, pp. 1918-1927 (2003), the contents of which are herein 
incorporated by reference). These aliasing artifacts are known as windmill artifacts. 
Presently, there are two types of methods to overcome this problem. One method changes 
hardware, and the other method alters software. 

[0005] The hardware modification applies a flying focus technique (in xy-plane, detector 
channel direction) to the z-axis direction and combines two projection sets with N samples in 
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z into one set of two N samples in z. (See T. Flohr, H.K. Bruder, K. Stierstorfer, S. Schaller, 
"Evaluation of Approaches to Reduce Spiral Artifacts in Multi-Slice Spiral CT 5 " RSNA 2003 
program, pp567, the contents of which are herein incorporated by reference.) However, this 
potential solution has the following disadvantages: (1) the data size becomes twice as large 
as currently used, (2) the scanners become more expensive, and (3) a new image 
reconstruction algorithm is required, among other things. 

[0006] The software solution may include non-adaptive z-filtering to reduce z-resolution 
(high frequency in z) or adaptive, object-dependent z-filtering. The non-adaptive solutions 
include z-filtering in the projection data domain with a fixed kernel (See K. Taguchi, U.S. 
Pat. No. 5,825,842 (1998), the contents of which are herein incorporated by reference), 
z-filtering in the image domain, and z-filtering in the projection data domain with a variable 
size of a kernel (See L Zmora, U.S. Pat. No. 6,560,308, the contents of which are herein 
incorporated by reference). 

[0007] One disadvantage of these solutions is that the z-resolution of the entire image is lost 
in a uniform fashion. Simply put, z-resolution is lost everywhere. Most of the causes of the 
windmill artifact, such as ribs, skulls, and spines, are located in the peripheral region of the 
images and the central region of images requires high spatial z-resolution. 
[0008] The adaptive z-filtering includes filtering in the projection data domain (See J. Hsieh, 
"Adaptive Interpolation Approach for Multi-slice Helical CT Reconstruction," Proc ofSPIE 
Vol. 5032, pp. 1876-1833 (2003), the contents of which are herein incorporated by reference) 
and filtering in the image domain. Both methods change the kernel of the z-filtering, 
depending on object dependent indices, such as the gradient in z and require a significant 
amount of processing, which is undesirable. 
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SUMMARY OF THE INVENTION 
[0009] In light of the above-described difficulties, the Applicant developed the present 
invention. A first non-limiting aspect of the invention provides a method for obtaining data 
from a computed tomography (CT) scan, including: obtaining projection data from at least 
two detector rows in a CT system; filtering the projection data in a direction of the at least 
two detector rows to obtain filtered data in which windmill artifacts are reduced; and 
reconstructing image data from the filtered data. 

[0010] Another aspect of the present invention provides an X-ray CT apparatus including: a 
helical scanning device configured to collect projection data while at least one of a gantry and 
a couch moves along an axial direction of the couch, the helical scanning device including, an 
X-ray source configured to generate X-rays, and a detector having detector elements arranged 
in at least two detector rows along the axial direction and configured to produce the 
projection data; and a processor including a filtering device configured to filter the projection 
data in a direction of the at least two detector rows to obtain filtered data in which windmill 
artifacts are reduced, and a reconstructing device configured to reconstruct the filtered data. 
[0011] Yet another aspect of the present invention provides an X-ray CT apparatus 
including: a helical scanning device configured to collect projection data while at least one of 
a gantry and a couch moves along an axial direction of the couch, the helical scanning device 
including, an X-ray source configured to generate X-rays, and a detector having detector 
elements arranged in at least two detector rows along the axial direction and configured to 
produce the projection data; and a processor including means for filtering the projection data 
in a direction of the at least two detector rows to obtain filtered data in which windmill 
artifacts are reduced, and a reconstructing device configured to reconstruct the filtered data. 
[0012] Another aspect of the present invention provides a computer program product storing 
instructions for execution on a computer system, which when executed by the computer 
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system, causes the computer system to perform the following steps: obtaining projection data 
from at least two detector rows in a CT system; filtering the projection data in a direction of 
the at least two detector rows to obtain filtered data in which windmill artifacts are reduced; 
and reconstructing the filtered data. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0013] A more complete appreciation of the invention and many of the attendant advantages 
thereof will be readily obtained as the same becomes better understood by reference to the 
following detailed description when considered in connection with the accompanying 
drawings, wherein: 

[0014] Figure la illustrates a multi-slice detector where N = 4; 
[0015] Figure lb represents geometry of the present invention; 
[0016] Figure 2a illustrates an example of a (r 2 D); 
[0017] Figure 2b illustrates a Gaussian distribution w(k)\ 
[0018] Figures 3a-3c illustrate the definition of a distance; 
[0019] Figure 4 illustrates an example of w Gn (0); 

[0020] Figure 5a illustrates an example of a function w z L used in both RF-soft and the RF- 
sharp algorithms; 

[0021] Figure 5b illustrates an example of a function wf used in the RF-soft algorithm; 
[0022] Figure 6 illustrates pixels projected onto the detector; and 
[0023] Figure 7 illustrates the z-resolution index. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0024] The sharpness of the kernel increases with decreasing r2D (e.g., a projected distance 
from the iso-center to the ray-sum onto the xy-plane), so that pixels near an iso-center have 



6 



• better z-resolution than in the periphery. Thus, the cause of the windmill artifact at the 
peripheral regions may be suppressed by applying a smoothing kernel to the corresponding 
detector channels. 

[0025] To this end, filter data in a detector row direction may be obtained prior to image 
reconstruction. The kernel of the filter may be defined as a function of the ray-angle. That is, 
the kernel may be adjusted based on the projected distance from the iso-center to the ray-sum 
onto the xy-plane. Equation 1 (below) provides the desired function: 

K 

Pout ( v > ch > row) = ^[w Cn (A:, ch) • p in (v, ch, row + k)] ( 1 ) 

k=-K 

In Equation l 9 p in is the projection data, v represents the projection number index 
corresponding to the projection angle /?, ch is the detector channel index corresponding to the 
ray angle 7, row is the detector row index corresponding to the cone angle a, w Gn is the 
coefficient of the z-filter, and 2K-1 represents the size of the kernel. 
[0026] Equations 2-4, below, represent the kernel: 

™ G n (*> ch) = w(k, ch) y ]T w(U ch) (2) 



-if i Y 

w(k,ch) = \/^xe _ K ^ K , (3) 



r 2D {ch) = R sin y(ch) = R sin{ l7m ' {Cch cK) 

\ Nch 



(4) 

J 



In Equations 2-4, R represents a distance from the focus to the iso-center, y m represents the 
maximum ray-angle, Nch represents the number of detector channels in one row, Cch is the 
detector channel index corresponding to 7= 0, v 2 d refers to the projected distance from the 
iso-center to the ray-sum onto the xy-plane, and o(r 2 D) defines a sharpness of the kernel. 
Figure 2a illustrates an example of <j(r 2D ), and Figure 2b illustrates w(k) using a different a. 
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[0027] Through this method, the sharpness of the kernel increases with a decreasing r^z), so 
that pixels near the iso-center achieve a better z-resolution than in the periphery. Thus, the 
windmill artifact at the peripheral regions is suppressed by applying the smoothing kernel to 
the corresponding detector channels (e.g., by smoothing filtering). 
[0028] The filtered projection data are then used for image reconstruction. Image 
reconstruction techniques include, but are not limited to a filtered backprojection and Fast 
Fourier Transform (FFT). According to the present invention, the 7-filtering may occur 
before the smoothing kernel is applied. 7-filtering may include applying a convolution 
process along a detector channel for image reconstruction such as ramp-filtering, but is not 
limited thereto. 

[0029] Figure 3 a shows the side view of a ray-sum, the rotation axis (z-axis), and the image 
to reconstruct (the voxel/pixel to backproject the datum/ray- sum). Figure 3b illustrates the 
two-dimensional distance, v 2 d, the projected distance from the iso-center to the ray- sum onto 
the xy-plane. Figure 3c illustrates the three-dimensional distance, r 3D , the distance from the 
iso-center to the voxel on the ray-sum where the ray-sum coincides with the xy-plane to 
reconstruct (particular z location). 

[0030] As a non-limiting alternative, Equation 5 (below) represents a case where the size of 
the kernel (2K+1) is fixed. 

™ Gn (*, ch)\ M = (1 - w Gn (0, ch)) l(2K) (5) 

An example of w Gn (Q) is represented in Figure 4. 

[0031] As a non-limiting alternative, it is possible to apply the filtering process in the 
detector row direction in the frequency domain rather than in the spatial domain. To this end, 
an FFT may be applied in the detector row direction, a channel-by-channel. The data having 
been transformed by the FFT may then be multiplied by a frequency representative of the z- 
filtering kernel. Finally, an inverse FFT may be applied to the filtered data. 
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[0032] Additionally, it is possible to apply a two-dimensional FFT by combining 7-filtering 
with row filtering according to the present invention. After the two-dimensional FFT is 
applied to the detector data, ramp-filtering may be applied in the 7-direction while smoothing 
filtering is applied in the row direction. 

[0033] As a non-limiting alternative embodiment, it is possible to filter data in a detector row 
direction prior to image reconstruction. The kernel of the filter may be changed, based on the 
distance from the iso-center to the pixel corresponding to the detector cell. The following 
Equation 6 represents the filtering equation: 



Pout ( V > Ck > WW) = Y,\~ W Gn (*> r 3D ) ' Pin ( V > Ck 9 VOW + k)] (6) 

k~-K 

[0034] Equations 7-16, shown below, describe the kernel: 

/ K 

w cn (k,r 3D ) = w(k, r 3D ) / £ w(i, r 3D ) (7) 



w(k, r 3D ) = \42n x e 



4— T 



■K<k<,K, 



(8) 



<r{r 3D ) = jr 3D /r 0 (9) 

r 3D =(r 3 D -r 0 )-wZ +r Q (10) 
' = {jx 2 +y 2 )= ^L 2 - 2/?z|cos/|+ R 2 (1 1) 

, ,v 2 Y m (Cch-ch) „„ N 

rich) = fm \ 1 (12) 

Nch 

L p ,a={L p , a -R)w L z +R (13) 



R + r„ 



R-r m < 
L 



<R + r„ 



R + r„ < 



<R~r„ 



(14) 
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L 



= zl tan a = 



Rz 



(15) 



d • (row — row c ) 



z = -CS-pi2n = -CS-(y-v c )/N^ 



vRev 



(16) 



[0035] In Equations 7-16, row c refers to the detector row index at G5=0, d is the projected 
detector cell height to the iso-center, v c denotes the projection number index at /3=0 (e.g., 
when the focus is in the plane to be reconstructed), CS is the table feed per rotation (helical 
pitch), N vRev is the number of projections per rotation. 

[0036] Two non-limiting variations are also possible when (3 is small: RF-sharp and RF-soft, 
where "RF" represents "Radius-Dependent Filtering Scheme." Using the RF-sharp 
algorithm, w z r is fixed at 1 .0. The w z L is a function of z(/3), as illustrated in Figure 5a. 
Equations 17 and 18 represent the RF-sharp algorithm. 

If w L 2 1, then L L(Y,p y a)\ r r(y 9 {3,a); w Gn -> w Gn (y,j3,a). (17) 

Else if w L z 0, then L R (fixed for all); r r(y)\ w Gn w Gn (y). (18) 

[0037] The RF-soft algorithm is a case where both w z L and w z T are a function of z(/3) as 
illustrated in Figures 5a-5b. Specifically, Figure 5a illustrates an example of a function w z L 
used in both the RF-soft and RF-sharp algorithms. Figure 5b illustrates an example of a 
function w/used in the RF-soft algorithm. Equations 19 and 20, below, represent the RF-soft 
algorithm. 



Else if w L z 0, then r -> r 0 (fixed for all); w Gn -> w Gn (y = y(r 0 ))( fixed for all). (20) 
[0038] Figure 6 illustrates the relationship between a detector cell and pixels to be 
reconstructed. The relationship illustrated in Figure 6 led the present inventors to the concept 
of r 3D . In Figure 6, the pixels were projected on the detector for /3=[-7r,7r]. The horizontal 
axis represents the detector channel, and the vertical axis represents the detector row. As 



If w r 2 -> 1, then r r(y,J3,a); w Gn -> w Gn (y,j3 9 a). 



(19) 
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illustrated in Figure 6, the projected region changes with j3. Figure 7 illustrates the z- 
resolution index using the RF-sharp algorithm and the RF-soft algorithm. 
[0039] All embodiments of the present invention conveniently may be implemented using a 
conventional general purpose computer or micro-processor programmed according to the 
teachings of the present invention, as will be apparent to those skilled in the computer art. 
Appropriate software may readily be prepared by programmers of ordinary skill based on the 
teachings of the present disclosure, as will be apparent to those skilled in the software art. 
[0040] As disclosed in cross-referenced U.S. Patent Application 6,236,05 1, a computer may 
implement the methods of the present invention, wherein the computer housing houses a 
motherboard which contains a CPU, memory (e.g., DRAM, ROM, EPROM, EEPROM, 
SRAM, SDRAM, and Flash RAM), and other optional special purpose logic devices (e.g., 
ASICS) or configurable logic devices (e.g., GAL and reprogrammable FPGA). The 
computer also includes plural input devices, (e.g., keyboard and mouse), and a display card 
for controlling a monitor. Additionally, the computer may include a floppy disk drive; other 
removable media devices (e.g. compact disc, tape, and removable magneto-optical media); 
and a hard disk or other fixed high density media drives, connected using an appropriate 
device bus (e.g., a SCSI bus, an Enhanced IDE bus, or an Ultra DMA bus). The computer 
may also include a compact disc reader, a compact disc reader/writer unit, or a compact disc 
jukebox, which may be connected to the same device bus or to another device bus. 
[0041] Examples of computer readable media associated with the present invention include 
compact discs, hard disks, floppy disks, tape, magneto-optical disks, PROMs (e.g., EPROM, 
EEPROM, Flash EPROM), DRAM, SRAM, SDRAM, etc. Stored on any one or on a 
combination of these computer readable media, the present invention includes software for 
controlling both the hardware of the computer and for enabling the computer to interact with 
a human user. Such software may include, but is not limited to, device drivers, operating 
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systems and user applications, such as development tools. Computer program products of the 
present invention include any computer readable medium which stores computer program 
instructions (e.g., computer code devices) which when executed by a computer causes the 
computer to perform the method of the present invention. The computer code devices of the 
present invention may be any interpretable or executable code mechanism, including but not 
limited to, scripts, interpreters, dynamic link libraries, Java classes, and complete executable 
programs. Moreover, parts of the processing of the present invention may be distributed 
(e.g., between (1) multiple CPUs or (2) at least one CPU and at least one configurable logic 
device) for better performance, reliability, and/or cost. For example, an outline or image may 
be selected on a first computer and sent to a second computer for remote diagnosis. 
[0042] [0004] The present invention may also be complemented with additional filtering 
techniques and tools to account for image contrast, degree of irregularity, texture features, 
etc. 

[0043] The invention may also be implemented by the preparation of application specific 
integrated circuits or by interconnecting an appropriate network of conventional component 
circuits, as will be readily apparent to those skilled in the art. 

[0044] The source of image data to the present invention may be any appropriate image 
acquisition device such as an X-ray machine, CT apparatus, and MRI apparatus. Further, the 
acquired data may be digitized if not already in digital form. Alternatively, the source of 
image data being obtained and processed may be a memory storing data produced by an 
image acquisition device, and the memory may be local or remote, in which case a data 
communication network, such as PACS (Picture Archiving Computer System), may be used 
to access the image data for processing according to the present invention. 
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[0045] Numerous modifications and variations of the present invention are possible in light 
of the above teachings. It is therefore to be understood that within the scope of the appended 
claims, the invention may be practiced otherwise than as specifically described herein. 
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